On the possibility of neutrino flavor identification at the highest energies 
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High energy astrophysical neutrinos carry relevant information about the origin and propagation 
of cosmic rays. They can be created as a by-product of the interactions of cosmic rays in the 
sources and during propagation of these high energy particles through the intergalactic medium. The 
determination of flavor composition in this high energy flux is important because it presents a unique 
chance to probe our understanding of neutrino flavor oscillations at gamma factors > 10 21 . In this 
work we develop a new statistical technique to study the flavor composition of the incident neutrino 
flux, which is based on the multipeak structure of the longitudinal profiles of very deep electron and 
tau neutrino horizontal air showers. Although these longitudinal profiles can be observed by means 
of fluorescence telescopes placed over the Earth's surface, orbital detectors are more suitable for 
neutrino observations owing to their much larger aperture. Therefore, we focus on the high energy 
region of the neutrino spectrum relevant for observations with orbital detectors like the planned 
JEM-EUSO telescope. 



I. INTRODUCTION 

The existence of neutrinos of E v > 10 12 eV is strongly 
motivated by the observation of cosmic rays up to en- 
ergies of order of 10 20 eV. These cosmic ray particles, 
mainly composed by protons and heavier nuclei, inter- 
act in the sources (see, for instance, Ref. [l|) or during 
their propagation through the intergalactic medium Q, 
producing weakly decaying particles such as pions and 
kaons that decay into neutrinos. High energy neutrinos 
can also be generated as the main product of the decay 
of superheavy relic particles 0, 0] • 

In general, the flavor content of the flux of astrophys- 
ical neutrinos is a function of the neutrino energy [5H7| . 
One of the main mechanisms that causes the flavor ra- 
tio to change is based on the fact that pions and muons 
produced in the sources lose their energy because of the 
interaction with the ambient magnetic field. The pion 
decay time is shorter than the corresponding one to the 
muon, and then, at sufficiently high energies, the prob- 
ability of pion decay, prior to significant energy loss, is 
larger than the corresponding one to the muon. As a 
consequence, the flavor ratio at the source changes from 
$u e ■ *„„ : $„ T = 1 : 2 : to : : -!',. 0:1: 0. 
The transition between theses two regimes depends on 
the characteristics of the source. Therefore, the deter- 
mination of the flavor content of the flux is of great im- 
portance for the understanding of the physical processes 
that take place in the source. The flavor information can 
also be used to study the more fundamental properties 
of neutrinos and to identify physics beyond the standard 
model (see, for instance, Refs. (^-IToj). 
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Cosmic ray observatories are sensitive to high energy 
neutrinos. In particular, there are two different ways 
to detect them. The first one consists in observing the 
development of horizontal air showers produced by the 
interactions of the neutrinos with the nucleons of the 
molecules of the Earth's atmosphere [ll|. The second 
one consists in observing the showers produced by the 
decay of taus generated by the interaction of tau neu- 
trinos that propagate through the interior of the Earth 

EMI. 

Horizontal and quasihorizontal neutrino showers as 
well as Earth skimming tau neutrino showers can be 
observed by using surface detectors or fluorescence tele- 
scopes placed over the Earth's surface or in space. At 
present no astrophysical neutrino has been observed in 
cosmic ray observatories and then, upper limits on its 
flux have been obtained [l5j |. 

In this work we develop a new technique to discrim- 
inate between scenarios with different compositions of 
electron and tau neutrinos, by using very deep tau and 
electron neutrino horizontal showers. This new technique 
is based on the morphology of the longitudinal profiles of 
the horizontal neutrino showers that can be observed by 
means of fluorescence telescopes. In particular, we study 
in detail the highest energy region of the neutrino flux 
that is relevant for space observations with the upcom- 
ing orbital fluorescence telescopes like JEM-EUSO [16} 
and the proposed Super- EUSO [ij} experiment. 



II. NEUTRINO SHOWERS 

Neutrinos can initiate atmospheric air showers when 
they interact with the nucleons of the air molecules. The 
probability that a neutrino interacts in the atmosphere 
increases with the zenith angle because of the increase of 
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the number of target nucleons. High energy neutrinos, 
propagating through the atmosphere, can suffer charge 
(CC) and neutral (NC) current interactions. The CC 
interactions are the most important for the space obser- 
vations because in the NC interactions most of the energy 
is taken by a secondary neutrino that could produce an 
observable air shower just in case it suffers a subsequent 
CC interaction. The observation of the showers produced 
by the hadronic component resultant from the NC inter- 
action depends on the energy threshold of the telescope. 
In this work just the showers initiated by CC interactions 
are considered. 

As a result of a CC interaction, a very high energy 
lepton, which takes most of the energy of the incident 
neutrino, is generated, 

n + N -> I + X, (1) 

where N is a nuclcon, I = {e~,/i~,r~}, and X is the 
hadronic component. Typically, the produced lepton I 
takes ~ 80% of the neutrino energy at E v = 10 20 eV, 
and the rest of the energy goes into the hadronic com- 
ponent X . In the case of electron neutrino showers, the 
produced electron together with the hadronic component 
initiate an air shower right after the neutrino interaction. 
Because most of the times the electron takes most of 
the neutrino energy, the electron neutrino showers have 
characteristics more similar to the electromagnetic show- 
ers than to the hadronic showers; however, the effects of 
the hadronic component are not negligible specially for 
showers developing in dense regions of the atmosphere 
(see Ref. 18] for details). 

In the case of tau neutrinos, while the hadronic com- 
ponent X initiates a low energy shower immediately af- 
ter the CC interaction (first bang), the tau lepton pro- 
duced propagates through the atmosphere almost with- 
out interacting and then, after a given distance, decays. 
The particles produced in the decay initiate a more ener- 
getic shower than the first one, producing a second bang. 
There are several channels for tau decay which can be 
classified in three groups flQj | . 

• Electromagnetic channel: t — > v T + e~ + v e - . 

• Hadronic channel: r — > v T + X, where X can be 
pions, kaons, etc. 

• Muonic channels: r — > v T + fjf + v^. 

The branching ratios of the electromagnetic, hadronic, 
and muonic channels are b ern = 0.18, = 0.645, and 
= 0.175, respectively. The tau neutrino showers cor- 
responding to the muonic channel are difficult to observe 
because of the large decay length of the muons at the 
energies considered, and then, only the electromagnetic 
and hadronic channels are relevant for the detection of 
tau showers. Therefore, considering just the electromag- 
netic and hadronic channels, ~ 78% of the tau neutrino 
showers are of hadronic origin whereas ^22% are of elec- 
tromagnetic origin. 



Muon neutrino showers are initiated by the electrons 
generated in the dominant decay channel of the muons 
produced in the CC interaction. However, these show- 
ers are difficult to observe because of the large decay 
length of the muon at such high energies, A(£^) = 
6.23 x 10 10 (i; A1 /10 19 eV) km. That is why muon neutrino 
showers are not considered in the subsequent analyses. 

Electron neutrino showers are simulated following the 
procedure introduced in Ref. [Hj], and the charge cur- 
rent neutrino-nucleon interaction is simulated by using 
PYTHIA 0, linked with the LHAPDF library [U, 
which allows the use of different sets of parton distribu- 
tion functions (PDFs). The CTEQ6 H3 set of PDFs, the 
most commonly used at the highest energies, is used in 
this work. Note that the use of a different set of PDFs has 
a small impact on the simulated neutrino showers fUjll^ . 
The secondary particles, produced in the interaction, are 
used as input in the program CONEX [24| (v2r2.3), in or- 
der to simulate the shower development. The high energy 
hadronic interaction model used for the shower simula- 
tions is QGSJET-II [25[. However, our main results are 
not strongly dependent on this assumption. 

The CC tau neutrino-nucleon interaction and the pro- 
duction of the corresponding tau lepton, are simulated 
using PYTHIA with the CETQ6 set of PDFs. The de- 
cay of the tau is simulated with the program TAUOLA 
[19J . Then, the particles produced in the decays are used 
as input in CONEX (with QGSJET-II) in order to sim- 
ulate the shower development. 

A library of horizontal neutrino showers is generated by 
using the simulation chain described above. The primary 
energy of the neutrinos ranges from 10 18 eV to 10 20 ' 5 eV 
in steps of A log(-E„/eV) = 0.25. The interaction point 
of the electron neutrinos and the decay of the taus are 
just above the Earth's surface, the densest region of the 
atmosphere. Note that the injection point is such that 
the trajectory of the showers starts at the vertical axis of 
a nadir-pointing orbital telescope at sea level. Therefore, 
hereafter we denote very deep horizontal (VDH) showers 
to the showers initiated by the interaction of horizontal 
electron neutrinos or the decay of tau leptons (generated 
by the interaction of horizontal tau neutrinos), such that 
the interaction point of the electron neutrinos, or the 
decay of the tau leptons, is placed at the vertical axis of 
a nadir-pointing telescope at sea level. 

Figure Q] shows the longitudinal profiles (deposited en- 
ergy per unit of atmospheric depth) for VDH electron and 
tau neutrino showers of E v = 10 20 eV. The multiple peak 
structure present in the profiles is caused by the fluctu- 
ations introduced by the Landau Pomeranchuk Migdal 
(LPM) effect, which is important for showers dominated 
by the electromagnetic component. Note that, as can be 
seen from the figure, the VDH electron neutrino showers 
are more affected by the LPM effect. 
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FIG. 1. Longitudinal profiles for VDH electron and tau neu- 
trino showers of E v = 10 20 eV. 



III. FLAVOR DISCRIMINATION 

As mentioned before, VDH neutrino showers present 
a multipeak structure caused by the fluctuations intro- 
duced by the LPM effect. In particular, electron neutrino 
showers are more affected by this effect because they are 
dominated by the electromagnetic component. In con- 
trast, tau showers are less affected by this effect, because 
~ 78% of the times the tau decays are of the hadronic 
type. Therefore, the number of peaks present in a given 
longitudinal profile should be a good parameter to dis- 
criminate between VDH electron and tau neutrino show- 
ers. 

Given a simulated longitudinal profile, the number of 
peaks and their positions are determined following the 
procedure described in Ref. [l8j |. The upper panel of 
figure shows the distribution of the position of the first 
peak, X^ nax , found from the beginning of the shower, Xq, 
for electron and tau neutrino showers of E v = 10 20 eV. 
Note that the distribution function of AJ^ al = X^ lax — 
Xq for VDH electron neutrino showers is bivalued and its 
first peak is located at AX^ al ~ 800 g cm -2 , while the 
second one is centered at AX^ ax ~ 1500 g cm -2 . The 
first peak is related to the development of the hadronic 
component of the cascade, whereas the second one mainly 



reflects the electromagnetic components of the shower 
(see Ref. [l8[ for a detailed discussion). As expected, 
the distribution function of AA^ oa . for VDH tau neutrino 
showers is also bivalued but the hadronic peak is much 
more important. The bottom panel of figure [2] shows the 
probability to find a profile with N{X* nax ) peaks. As 
expected, the profiles corresponding to electron neutrino 
showers have a larger probability to have more than one 
peak. In particular, the probability to have just one peak 
for VDH electron neutrino showers of E v = 10 20 eV is ~ 
0.67, whereas for VDH tau neutrino showers it is ~ 0.98. 
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FIG. 2. Upper panel: distribution of the first peak AXmax 
for VDH electron and tau neutrino showers. Bottom panel: 
probability to find N(Xl nax ) peaks in a profile corresponding 
to a VDH electron and tau neutrino showers. The primary 
energy of the neutrinos is E v = 10 20 eV. 

Figure |3] shows the probability to find just one peak, 
Px 1 i m a profile corresponding to VDH electron and 
tau neutrino showers as a function of the primary energy. 
The solid lines in the figure are fits to the data points with 
the function 

Po 
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where po,E , and 7 are free fit parameters. The proba- 
bility to find a profile with just one peak decreases with 
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primary energy, owing to the fact that the LPM effect be- 
comes more important at higher energies. As expected, 
Px~l (E v ) for tau neutrinos is larger than the corre- 
sponding one for electron neutrinos. 




FIG. 3. Probability to find a profile with just one peak as a 
function of primary energy for VDH electron and tau neutrino 
showers. 
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The decay length of a tau lepton of energy E T is given 



L T {E T ) = 4900 



10 20 eV 



km. 



(3) 



Therefore, the probability that a high energy tau, gener- 
ated in a CC interaction of a horizontal tau neutrino with 
an atmospheric nucleon, decays just above the Earth's 
surface and is smaller than the probability that an hori- 
zontal electron neutrino suffers a CC interaction with an 
atmospheric nucleon at the same point. Because of this, 
the number of VDH tau neutrino showers are suppressed 
compared with the electron neutrino ones at the higher 
energies. In particular, the ratio between the number of 
VDH tau neutrino showers and VDH electron neutrino 
showers for an incident flux with the same number of 
electron and tau neutrinos is given by 

dP sh I dP sh 

Ru(Eu) = -g-{E v ,t = Q)j-^-{E v ^ = 0) 

= expWW )) r dErP{EAEv) 

Po Jo 



x / deexp(-*(0/A„(^)) p(h(0) 
Jo 

exp(-t/L T (E r )) 



Lr(E T ) 



(4) 



Here £ is the distance in the horizontal direction mea- 
sured from a point just above the Earth's surface, 
dPff I 'd^{E v X = 0) is the probability per unit of length 
to initiate an horizontal /-flavor neutrino shower at £ = 0, 
P(E T \E V ) is the probability to generate a tau lepton 
of energy E T given the CC interaction of a tau neu- 
trino of energy E v with a nucleon of the atmosphere, 



\ V (E V ) = 7.23 x 10 7 (^/I0 19 eV)-°- 358 g cm- 2 is the 
mean free path for CC neutrino-nucleon interaction cal- 
culated using the cross section of Ref. [2a |. p(h) is the 
atmospheric density as a function of the altitude h, 



is the altitude as a function of £ with i?g 
the Earth, 



x(0 



(5) 

the radius of 



(G) 



and X = X(£ = 0). 

Figure [4] shows R v as a function of the neutrino en- 
ergy. The atmospheric density profile used in the calcu- 
lation corresponds to Linsley's parametrization of the US 
standard model l 27], and the probability P(E T \E V ) is ob- 
tained from PYTHIA simulations. R v decreases with the 
neutrino energy in such a way that between E v = I0 18 
eV and E v = 10 19 75 eV it falls from ~ 1 to ~ 0.13. 
Note that when the tau neutrino energy decreases, the 
energy of the tau lepton produced in the interaction also 
decreases; as a consequence the lifetime of the tau lep- 
ton, in the laboratory frame, also decreases. Therefore, 
when the tau neutrino energy is small enough, the tau 
lepton decays almost immediately after being produced, 
and then the shower starts at £ = as for the case of 
electron neutrinos, producing values of R v close to one. 
In any case, for L T — > the function exp(— £/L r )/L r , in 
Eq. (fJJ, can be approximated by the Dirac delta func- 
tion, <5(£). Then it can easily be seen that i?„ = 1 when 
the Dirac function is used in Eq. (QJ. 
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FIG. 4. Ratio of the number of VDH tau neutrino showers to 
the number of VDH electron neutrino ones as a function of 
primary energy. 

Although the number of VDH tau neutrino showers is 
smaller compared with the VDH electron neutrino ones, 
few events can give relevant information about the dif- 
ferent flavors present in a sample. Figure [5] shows the 
regions of 95% of probability to find a fraction of neu- 
trino showers with one peak, fx — ni/N, as a func- 
tion of the sample size N for VDH electron and tau 
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neutrino showers of E v — 10 1975 eV. For the case of 
VDH electron neutrino showers this region is enclosed by 
fx = and /i = f{ nax (N) and for tau neutrino show- 
ers by fx = f™ m (N) and fx = 1 where f? ax (N), and 
f™ m (N) are calculated by using nx as a binomial ran- 
dom variable. The probability to find just one peak in 
a longitudinal profile is obtained from the fits of P^i 
(see Eq. ©). It can be seen that with a sample of a small 
number of events, it is possible to say something about 
the neutrino flavors present in the sample, depending on 
the observed number of profiles with just one peak, nx- 
If fx, obtained from the observations, falls in the non- 
shadowed region, the hypothesis of having a pure sample 
of electron or tau neutrinos is rejected with at least 95% 
probability. Also for N > 10, if fx falls in some of the 
two shadowed regions, it indicates that it is possible to 
reject the hypothesis of having a pure sample of neutri- 
nos of the opposite flavor, with at least 95% probability. 




Here the factor 1 — takes into account the fact that the 
tau neutrino showers that come from the muonic channel 
of the tau decay are not considered. For a given value 
of the sample size N, the number of VDH electron neu- 
trino showers is obtained by taking at random a number 
from the binomial distribution of Eq. ([7]). If in the ith 
step of the simulation the number of VDH electron neu- 
trino showers is n\, and the number of VDH tau neutrino 
showers is n\,^ = N — , the number of profiles with 
just one peak is obtained by sampling the binomial dis- 
tributions B{n\, ,p\, ) and B{n l v ^,p x VT ) where pi and p\ 
are obtained from the fits of Px 1 as explained before. 

Figure [B] shows the intervals of 68% of probability of fx, 
as a function of the sample size, for three different values 
of electron neutrino abundances, c„ c = {0.25,0.5,0.75}, 
and for E v — 10 19 - 75 eV. It can be seen that scenarios 
in which the incident flux is dominated by tau neutrinos 
are easier to discriminate, i.e. samples of a smaller size 
are required. Given the flavor content of the incident 
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FIG. 5. Regions of 95% probability of the random variable 
fi = ni/N (see the text) as a function the sample size N for 
VDH electron (red) and tau (green) neutrino showers with 



The distribution function of fx for a given energy, sam- 
ple size, and the electron neutrino abundance of the in- 
cident flux for a binary mixture of electron and tau neu- 
trinos, c Ve — N(y e )/(N(v e ) + N(v T )), is obtained from 
Monte Carlo simulations. Assuming that the distribu- 
tion function of the number of VDH neutrino showers of 
a given flavor I, n Vl , corresponding to a given observa- 
tion time is Poissonian, the probability to find n„ c and 
n UT showers in a sample of size N — n v<s + n Vr is given 
by the binomial distribution 

P sh {n Ue \N) = B(N,p) = V"«(l -p) N - n »°, (7) 
where 

P(Eu) = — c • (8) 

1 + — (l-M Rv{E v ) 



FIG. 6. Regions of 68% probability of /i for three different 
values of electron neutrino abundance, c„ e — {0.25, 0.5, 0.75}. 
The primary energy corresponds to E v — if) 19,75 eV. 

flux, it is possible to use the parameter fx to calculate 
the minimum number of events, N m i n , required to reject 
a false hypothesis. The notation used for the different 
scenarios considered is $ t/e : $„ T = 1 : 1 for the case 
in which there are an equal number of electron and tau 
neutrinos on the incident flux, $ I/e : $> Vt =0:1 means 
that there are just tau neutrinos, and §> Ue : <& Ut = 1:0 
means that there are just electron neutrinos. Figure [7] 
shows an example of the calculation of N m i n for the case 
in which there are just tau neutrinos on the incident flux 
and for E v = 10 19 75 eV. The shadowed region corre- 
sponds to the region of 68% probability of fx for the case 
in which there are just tau neutrinos in the sample, i.e., 
0:1. The constant /i = together with the black line 
and also fx — together with the red solid line, enclose 
the regions of 95% probability to find fx for the cases 
in which there are an equal number of tau and electron 
neutrinos, 1:1, and just electron neutrinos, 1 : 0, on the 
incident flux, respectively. Therefore, in the 68% of the 
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cases it is possible to reject the hypothesis 1 : with 95% 
probability with a sample of size N > N^ n = 31 and 
also, in the 68% of the cases, it is possible to reject the 
hypothesis 1 : 1 with 95% probability by using a sample 
of size N > N m : } n = 33. Note that, if the region of 95% 
probability of /i, corresponding to : 1, is used instead 
of the region of 68% probability, larger values of N m in 
are obtained. 




FIG. 7. /i as a function of N for E v = 10 19 ' 75 eV. Shadowed 
region shows the 68% probability of the distribution function 
of /i for the incident flux with flavor content 0:1. N^} n and 
N^n correspond to the minimum number of events required 
to reject the hypothesis 1 : 1 and : 1 with 95% probability 
(see text for details). 

Figure [8] shows N m in as a function of primary energy 
for three different hypotheses of the flavor content of the 
incident flux. The rejection probability is 95% and the 
calculation is done considering the regions of 68% prob- 
ability (blue squares) and 95% probability (red circles) 
of the distribution function of f\ for the true hypothesis. 
N m in depends on the difference between the probabilities 
to find a profile with just one peak corresponding to VDH 
electron and tau neutrino showers, Api = p 1 ^ — p 1 ^, and 
also on the parameter i?„. The larger the values of Api 
and the closer to one the values of R v the easier to dis- 
criminate among different scenarios. From figures [3] and 
2] it can be seen that while Api increases with the neu- 
trino energy, R v decreases because of the increase of the 
tau decay length with primary energy. As a consequence, 
samples with larger size are required. 

The top panels of figure [5] show the case in which the 
true hypothesis is 1 : 1. It can be seen that the val- 
ues of N m i„ required to reject the hypothesis of having 
tau neutrinos alone decreases with primary energy. This 
happens mainly because Api increases with the neutrino 
energy and also, but less important, because the distri- 
bution function of fx, for the 1 : 1 case, becomes more 
similar to the 1 : case for increasing values of the pri- 
mary energy, which is caused by the decrease of ¥L V with 
E v . To test the hypothesis of having electron neutrinos 
alone is more difficult for increasing values of E v , be- 



cause, as said before, the distribution function of /i for 
the 1 : 1 case becomes more similar to the corresponding 
one to the case in which there are electron neutrinos alone 
on the incident flux. This effect is more important than 
the increase of Api with E v causing N m in to increase 
with the neutrino energy. Note that the values of N m in 
required to reject the hypothesis of having tau neutrinos 
alone is more than 2 orders of magnitude smaller than 
the ones required to reject the hypothesis of having elec- 
tron neutrinos alone. In particular, for E u — 10 19 7 eV 
and considering the regions of 95% probability of f\ for 
the true hypothesis, N min = 45 to reject the hypothesis 

: 1 with 95% probability and N min = 5529 to reject the 
hypothesis 1 : with 95% probability. 

The middle panels of figure [5] show the case in which 
the incident flux is composed by electron neutrinos only, 

1 : 0. For similar reasons discussed for the previous case, 
Nmin is a decreasing function of energy when the hypoth- 
esis : 1 is considered. N m in is an increasing function of 
energy when the hypothesis 1 : 1 is considered and N m in 
is more than 2 order of magnitude larger compared with 
the corresponding one to the test of hypothesis 0:1. 

Finally, the bottom panels of figure [5] show the scenario 
in which the incident flux is composed by tau neutrinos 
alone. In both cases N m in is a decreasing function of 
energy because of the increase of Api with primary en- 
ergy. As expected, to reject the hypothesis 1 : 1 requires 
a larger number of events because of the contribution of 
tau neutrinos to f\. In any case, at E u = 10 18 eV for 
95% of the cases, it is possible to reject both hypotheses 
with < 200 events and with rejection probability of 95% 
events. At E v = 10 20 eV just < 20 events are required. 

Note that the energy deposition of electron and tau 
neutrino showers is very different because the tau lep- 
tons, generated in the CC interactions of tau neutrinos 
with the nucleons of the atmosphere, decay to tau neu- 
trinos and other particles; thus the energy taken by this 
secondary tau neutrino do not go to the shower. There- 
fore, if the deposited energy is used to reconstruct the pri- 
mary energy without any correction, the standard proce- 
dure when fluorescence telescopes are used to observe the 
showers, then tau showers are reconstructed with smaller 
energies than that corresponding to electron neutrinos. 
The impact of the energy determination on the method 
developed above is studied in detail in Appendix. 



IV. CONCLUSIONS 

In this work we have developed a new statistical tech- 
nique intended to study the flavor content of the inci- 
dent flux of high energy astrophysical neutrinos. This 
new method is based on the morphological differences 
between the longitudinal profiles corresponding to very 
deep horizontal electron and tau neutrino showers that 
can be observed by means of fluorescence telescopes. 
In particular, the multipeak structure of the profiles 
strongly depends on the flavor of the incident neutrino. 
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FIG. 8. Minimum number of events, N m in, as a function of primary energy, needed to reject different hypotheses given true 
scenarios. The red circles correspond to the region of 95% probability of /i for the true hypothesis and the blue squares to the 
region of 68% probability. The rejection probability considered is 95%. 



Very deep horizontal showers initiated by electron neutri- 
nos present more fluctuations, owing to the LPM effect, 
than the corresponding one to tau neutrinos. Therefore, 
the number of showers with just one peak can be used 
to discriminate between different scenarios for the flavor 
content of the incident neutrino flux. 

The difference between the probability to find just one 
peak, in a given profile, for very deep horizontal electron 
and tau neutrino showers starts to be important for ener- 
gies > 10 19 ' 5 eV. Then this technique is more relevant for 



fluorescence telescopes in orbit around the Earth, like the 
upcoming JEM-EUSO mission and Supcr-EUSO, owing 
to their huge exposure. 

Scenarios in which the incident flux is dominated by 
tau neutrinos are easy to identify, and samples with sev- 
eral events, depending on the primary energy, are re- 
quired in order to reject the hypothesis of having electron 
neutrinos alone. However, because of the suppression of 
very deep horizontal tau neutrino showers with respect 
to electron neutrino ones, a large number of events is 
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required to discriminate between scenarios with electron 
neutrinos alone and a mix of an equal number of tau and 
electron neutrinos. 
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Appendix A: Energy deposition of neutrino showers 

Tau neutrino showers are initiated by the decay of the 
tau leptons, generated in CC interactions of tau neutri- 
nos with atmospheric nucleons. As mentioned before, the 
generated tau lepton takes part of the energy of the inci- 
dent tau neutrino. Also, in all decay modes of the taus, 
a tau neutrino is generated that takes part of the energy 
of the tau. This energy is missed because the neutrinos 
do not participate in the shower development. Figure [§] 
shows the median and the region of 68% probability of 
the deposited energy as a function of the incident neu- 
trino energy for VDH tau and electron neutrino showers. 
The deposited energy is obtained by integrating the lon- 
gitudinal profiles of the showers. The deposited energy 
is the observable used to infer the energy of the primary 
particle when the fluorescence telescopes are used to ob- 
serve the air showers. From the figure it can be seen that 
energy deposition of tau neutrino showers is smaller than 
the corresponding one to electron neutrino showers with 
the same primary energy. Also the fluctuations are much 
larger for the case of tau neutrinos. 
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FIG. 9. Deposited energy of VDH electron and tau neutrino 
showers as a function of primary energy. The solid lines are 
the median of the distributions and the shadowed regions cor- 
respond to the regions of 68% probability. 



The top panel of figure [10] shows the distribution of the 
ratio between the deposited energy, Ejj , and the tau neu- 
trino energy, E Ut , for E Vt = 10 18 eV and E Vt = 10 20 5 
eV, the minimum and the maximum values of energy con- 
sidered, respectively. It can be seen that the shapes of 
these distributions do not change significantly with the 
tau neutrino energy. The large fluctuations on the de- 
posited energy of tau neutrino showers are dominated by 
the fluctuations of the energy taken by the tau neutrino 
generated in the tau decays. The latter can be inferred 
from the plot at the bottom panel of figure [TUJ which 
shows the distribution of the ratio between the energy 
that goes to the shower, E s h, and the tau neutrino en- 
ergy for the same values of E Vt considered before. E s h 
is obtained by adding the energy of the particles gener- 
ated in the tau decays except the one of neutrinos. It 
can be seen that the energy that goes to the showers 
presents large fluctuations and also that, on average, the 
deposited energy is smaller than E s h, 8% for E Vt = 10 18 
eV and 2% for E Vr = 10 20 5 eV. This is due to the fact 
that there is also missing energy in the development of 
the shower. 

For the case of electron neutrinos the deposited energy 
is smaller than the primary energy by less than ~ 1%, 
and the fluctuations, <r[£ , £)]/med(S£)), are smaller than 
~ 0.3% in the whole energy range under consideration. 

As mentioned before, the deposited energy of the show- 
ers is used to infer the primary energy. Thus, without 
any correction, the reconstructed energy of tau neutrino 
showers is smaller than the corresponding one to electron 
neutrino showers. On average, the deposited energy of 
tau neutrino showers goes from ~ 0.4 E Ut at E Vt = 10 18 
eV to - 0.5 E Vt at E Vt = 10 20 5 eV. Therefore/because 
of the fact that P^ 1 decreases with the neutrino en- 
ergy, the difference between Px 1 for tau and electron 
neutrinos becomes smaller when the deposited energy is 
used to estimate the neutrino energy, decreasing the dis- 
crimination power of f\. 

The distribution of Ed/E Ut can be fitted by a third 
degree polynomial, as can be seen from the top panel of 
figure [TDJ Because of the small variation of the shape of 
this distribution with E„ T , the average between the fits 
of the distributions corresponding to E Vr = 10 18 eV and 
E Ut = 10 20 ' 5 eV is used for the whole energy range under 
consideration. The solid line on the plot of the top panel 
of figure [TUl corresponds to the average function, denoted 
as P(E d \E Ut ). 

The probability to find just one peak in a given profile 
as a function of the deposited energy is given by, 



P xlm jE D )= dE Vr P XL {E Vr )P{E Vr \E D ) (Al) 
Jo 



where 



P{E Vr \E D ) 



P(E D \E Ur ) P{E Vr ) 
Io° dE UT P(En \E Vt ) P(E Vt ) 



(A2) 



Here, the energy distribution of the neutrinos is taken as 
a power law, P[E Vt ) = KE~i . 
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log(E /eV) = 20.5 




difference with the ideal case. 



FIG. 10. Top panel: distribution of the ratio between the 
deposited energy and the tau neutrino energy for E Vt — 10 18 
eV and E„ T — 10 20 ' 5 eV. Bottom panel: distribution of the 
ratio between the energy that goes to the shower (see text) 
and the tau neutrino energy for E Vt = 10 18 eV and E Vt = 



eV. Dashed lines are fits of the distributions with third 



degree polynomials, and the solid line in the top panel is the 
average of the two fits. 



For the case of electron neutrinos, the approximation 
P(E D \E Ue ) = S(E D - 0.99 E Ve ) is used. Here 5(x) is the 
Dirac delta function. By using Eq. (|A1[) it is easy to 
show that P x i ma!c (E D ) ~ Px^JEd /0.99), 

Figure Qj] shows Px} na7 . (Ed) for VDH electron and tau 
neutrino showers as a function of the deposited energy. 
The solid lines correspond to the ideal case in which 
Ed = E v ; i.e., they correspond to the fits shown in fig- 
ure [3] It can be seen that, for VDH electron neutrino 
showers, P x i {Ed) is practically indistinguishable from 
P x i {Ed)- However, for the case of VDH tau neutrino 
showers Px 1 {Ed) is smaller than the ideal case and the 
difference increases with the deposited energy. Note that, 
for the case of VDH tau neutrino showers, Px^ (Ed) 
depends on the energy distribution of the incident flux 
(see Eq. (|A2[> ). From the figure it can be seen that, as 
expected, the steeper the neutrino flux the smaller the 




FIG. 11. P x i as a function of log(Su/eV). The red lines 
correspond to VDH electron neutrino showers and the green 
lines to VDH tau neutrino showers. The solid lines correspond 
to the ideal case in which Ed = E v . The dashed red line 
corresponds to VDH electron neutrino showers for which the 
approximation En = 0.99 E v is used. 

The ratio between the number of VDH tau neutrino 
showers and the number of VDH electron neutrino ones, 
as a function of the deposited energy, is given by (see Eq. 
@) 



Rv(E v ) 



dP. 3 . 



-(EdA = 0) 



'dP, 3 . 



d£ 



-(E D ,£ = 0) (A3) 



where 

d psh 



ish poo dP s 

^(E D ,0) = / dE„ T -^(E UT) 0)P(E UT \E D ),(A4) 
£ Jo <% 

(A5) 



dt, jo 

dPf, dP» h , 

-^(E D ,0) = -£-(E D /0M,0) 



Figure [T2l shows R v (Ed) as a function of the deposited 
energy for 7 = 2 and 7 = 3. For log(-Eo/eV) > 18.5, 
R V (E) is smaller than R V (E). 

Figure [T3] shows JV mm as a function of the deposited 
energy for 7 = 2 obtained from simulations following the 
procedure describe in Sec. Mil 

From figure [151 it can be seen that N min (E) is larger 
than N mm (E). The ratio r(E) = N min (E)/N min (E) can 
be approximated by a linear function of the logarithm of 
the deposited energy. Table Q] shows the variation of r 
from E D = 10 19 2 eV to E D = 10 20 5 eV obtained form 
the linear fits for the different assumptions about the in- 
cident flux and for the corresponding tests. The results 
shown in Table U are obtained considering the region of 
95% probability for the true hypothesis, and similar re- 
sults are obtained when the region of 68% probability is 
considered. It can be seen that the discrimination be- 
tween 1 : 1 and 1 : is more affected by the use of the 
deposited energy instead of the true energy, which is due 
to the smaller values obtained for R y (E) compared with 
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FIG. 12. Ratio between the number of VDH tau neutrino 
showers and the number of VDH electron neutrino ones as a 
function of the deposited energy. 



the corresponding ones for i? t/ (_E). The discrimination 
between these two scenarios is very sensitive to i?„ and 
to R u when the deposited energy is considered. 



TABLE I. r(10 192 eV) -> r(10 20 5 eV) for 7 = 2 and 7 = 3, 
obtained considering the region of 95% probability for the 
true hypothesis. 



True Test 7 = 2 7 = 3 

T7I TTo 3.1 -> 3.2 1.8 ->• 2.4 

1:1 0:1 1.8->1.2 1.2 1.2 

1:0 1:1 3.2^3.2 1.8 2.4 

1:0 0:1 1.8^1.4 1.2 -»■ 1.4 

0:1 1:1 1.9 1.0 1.3 1.2 

0:1 1:0 2.0^1.0 1.4 -»■ 1.0 
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FIG. 13. Minimum number of events, N m in, as a function of the deposited energy, needed to reject different hypotheses given 
true scenarios. The energy distribution of the incident neutrino flux corresponds to a power law of spectral index 7 = 2. 



